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Piezo-magnetoelectric effect, namely simultaneous induction of both the ferromagnetic moment 
and electric polarization by an application of uniaxial stress, was achieved in the non-ferroelectric 
and antiferromagnetic ground state of DyFeOs. The induced electric polarization and ferromagnetic 
moment are coupled with each other, and monotonically increase with increasing uniaxial stress. 

The present work provides a new way to design spin-driven multiferroic states, that is, magnetic 
symmetry breaking forced by external uniaxial stress. 

PACS numbers: 75.85.+t,74.62.Fj 


The discovery of the spin-driven ferroelectricity in 
TbMnOsJH opened the door to new frontiers in multi- 
ferroics, in which different kinds of ferroic orders coexist 
@0. The extensive studies on recently-explored multi- 
ferroics, such as orthorhombic RMnOs (R is a rare-earth 
element) [l|,[B, 6 ], M 3 V 2 O 8 0 , MnWCU Q etc., demon¬ 
strated that ferroelectricity can arise from the inver¬ 
sion symmetry breaking due to magnetic orders. As for 
the microscopic mechanisms to explain the spin-driven 
ferroelectricity, the most robust scheme is the inverse 
Dzyaloshinskii-Moriya (IDM) mechanism [9], in which 
the magnetically induced electric polarization (P) is de¬ 
scribed as P oc eij x (Si x Sj ) [10, EH , where is 
a unit vector connecting between two spins, Si and Sj. 
This formula predicts ferroelectricity in a ‘cycloidal’ spin 
order, in which the spontaneous electric polarization ap¬ 
pears along the twofold axis perpendicular to both the 
propagation wave vector and (Si x Sj ), regardless of the 
symmetry of the underlying chemical lattice. 

On the other hand, there are two other mechanisms for 
the spin-driven ferroelectricity 12-151, namely, the spin- 
dependent p-d hybridization model |l3| and the magne¬ 
tostriction model 0, 0. In these models, the emer¬ 
gence of the ferroelectricity cannot be explained by local 
spin arrangement alone, but is explained by taking into 
account the symmetry of the crystal structure. This sug¬ 
gests that in the non-IDM type spin-driven multiferroics, 
the ferroelectricity can be controlled via changes in crys¬ 
tal structural symmetry. 

One of the most primitive ways to control the symme¬ 
try of the crystal structure is an application of anisotropic 
stress. It has recently been demonstrated that an ap¬ 
plication of compressive uniaxial stress can induce spin- 
driven electric polarization in Ba 2 CoGe 207 0 , in which 
the origin of the ferroelectricity is explained by the p-d 
hybridization model 17]. However, this system originally 
has a non-centrosymmetric (but non-polar) crystal struc¬ 
ture belonging to the space group of P42ira, in which the 
application of uniaxial stress naturally leads to piezoelec¬ 
tric polarization regardless of the spin degree of freedom. 
Therefore, the uniaxial-stress-induced ferroelectricity in 


the antiferromagnetic phase might be trivial. 

In the present study, we have investigated uniaxial- 
stress effects on magnetic and dielectric properties of 
a dysprosium ortho ferrite, DyFeOs, which has a cen- 
trosymmetric crystal structure. This system is known to 
display multiferroic nature in a magnetic-field-induced 
phase above ~ 3 T [II|. The multiferroicity can be 
described in terms of the magnetostriction between the 
magnetic moments on Fe and Dy atoms 00. On 
the other hand, the ground state is an antiferromagnetic 
and non-ferroelectric state. In this Letter, we report 
that an application of uniaxial stress induces ferroelec¬ 
tricity in the ground state. Moreover, the uniaxial-stress- 
induced ferroelectricity is accompanied by (weak) ferro¬ 
magnetism. To the best of our knowledge, this is the 
first experimental observation of ‘piezo-magnetoelectric 
effect’, namely simultaneous induction of both the elec¬ 
tric polarization and ferromagnetic moment, in a single¬ 
phase compound. The present study also reveals that the 
induced ferromagnetic moment and electric polarization 
are tightly coupled with each other, indicating that the 
piezo-magnetoelectric responses are of spin origin. 

DyFeOs has an orthorhombically-distorted perovskite 
structure, as shown in Fig. QJb). In this Letter, we 
employed a conventional Pbnm (orthorhombic) setting. 
The magnetic property of this system has been investi¬ 
gated since 1960s 0 l2f| . At room temperature, the 
system exhibits the most conventional-type (G-type) an- 
tiferromagnetic order, in which the staggered Fe 3+ spins 
direct nearly along the a axis, while the magnetic mo¬ 
ments of Dy 3+ ions are disordered. Because of the DM 
interaction between the Fe spins, the spins are slightly 
canted from the a axis, yielding a parasitic layered (A- 
type) antiferromagnetic and weak ferromagnetic compo¬ 
nents along the b and c axes, respectively [2l| . According 
to the Bertaut’s notation |22|, this magnetic structure is 
described as G x A y F z . We refer to this phase as the weak- 
ferromagnetic (WFM) phase. 

With decreasing temperature in zero magnetic field, 
the system undergoes a magnetic phase transition from 
the WFM phase to an antiferromagnetic (AFM) phase 
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(b) Dy-ordered AFM phase, [F e{A x ,G y ,C z ), Dy (G x ,A y )] 
Fe (z=0), Dy (z=0.25) Fe (z=0.5), Dy (z=0.75) 
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FIG. 1: (Color online) (a) Crystal structure of DyFeOs. Pur¬ 
ple dotted lines show the paths of the nearest neighbor ex¬ 
change coupling between the Fe spins and Dy moments. Mag¬ 
netic structures of (b) the Dy-ordered AFM phase and (c) the 
field-induced multiferroic phase. 


around Xsr = 40 K. This transition was identified to be 
a reorientation of the majority G-component from along 
the a axis to the b axis [20|. This spin-reorientation leads 
to disappearance of the WFM component, and instead, 
the canted spin components, namely the a and c compo¬ 
nents exhibit the A-type and G-type orders, respectively. 
As a result, the magnetic order is described as A x G y C z . 
This indicates that the direction of the G-component is 
relevant to the weak ferromagnetism in this system. 

With further decreasing temperature, Dy moments 
also exhibit an antiferromagnetic order below T^ y = 4 
K. We refer to this phase as the Dy-ordered AFM phase. 
According to a previous magnetization study on isostruc- 
tural DyAlC >3 [23j , the configuration of the Dy moments 
was reported to be G x A y . In Fig. DJc), we show the spin 
arrangements in the Dy-ordered AFM phase. When a 
magnetic field is applied along the c axis in this phase, 
the system exhibit a magnetic-field-induced phase transi¬ 
tion, at which the configuration of the Fe spins turns back 
to G x A y F z , around H$r ~ 3 T. Tokunaga et al have 
reported that the electric polarization appears along the 
c axis in the field-induced phase [1§|]. They explained the 
origin of the multiferroic nature by adapting the magne¬ 
tostriction model to the nearest neighbor exchange cou¬ 
pling between the Fe spins and Dy moments, as shown 

in Fig. He) [ 13 , 0 . 

Piezoelectric and piezomagnetic responses for these 
magnetic phases can be predicted in terms of the mag¬ 
netic point group theory jioj. For example, the mag¬ 
netic point group for the AFM phase and the Dy-ordered 
AFM phase are ‘mnmnrrC and ‘222’, respectively. By apply¬ 
ing uniaxial stress along the (110) direction in the AFM 


phase, we can break two mirrors out of the three, and 
therefore the magnetic point group is reduced to c 2/m\ 
which arrows ferromagnetic moment to appear along the 
c axis. Similarly, in the Dy-ordered AFM phase, an appli¬ 
cation of uniaxial stress along the (110) direction reduces 
the magnetic point group symmetry down to ‘2’, lead¬ 
ing to ferroelectric polarization as well as ferromagnetic 
moment parallel to the c axis. To verify these symmetry 
arguments, we performed magnetization and pyroelectric 
measurements under applied uniaxial stress. 

A single crystal of DyFeOs was grown by the floating- 
zone method and cut into a rectangular shape with di¬ 
mensions of 2.1 x 2.3 x 0.9 mm 3 . The widest surfaces 
were normal to the c axis, and two of the other sur¬ 
faces were selected to be the (110) plane. The magne¬ 
tization measurements were performed using Magnetic 
Property Measurement System (Quantum Design Inc.) 
with the uniaxial-stress insert used in the previous stud¬ 
ies H3 0. For this experiment, we also developed 
a horizontal-uniaxial-stress cell, which is essentially the 
same as that used in the previous pyroelectric measure¬ 
ments on Ba 2 CoGe 207 [16|. The uniaxial compressive 
stress cr was applied on the (110) surfaces. Magnetic 
fields H were applied along the c axis. For the pyroelec¬ 
tric measurements, silver-paste electrodes were formed 
on the surfaces normal to the c axis in order to observe 
the electric polarization along the c axis (P c ). We also 
employed the uniaxial-stress insert for Physical Property 
Measurement System (Quantum Design Inc.), which was 
used in the previous studies |16|, |24j, [25| . We measured 
displacement electric current with varying temperature 
T or a using an electrometer (Keithley 6517B). By inte¬ 
grating the current with respect to time, we deduced P c . 
Note that the magnitudes of cr were calculated from the 
load applied from the top of the insert and the pressur¬ 
ized area of the sample. 

Figure [2fa) shows the results of the magnetization 
measurements at 6 K in the AFM phase. We observed 
that the M-H curve shows a distinct hysteresis around 
H = 0 under cr, revealing that the application of cr in¬ 
duces weak ferromagnetism in the AFM phase. On the 
other hand, the critical field between the AFM phase to 
the field-induced WFM phase is hardly affected by the 
application of cr. 

We also measured magnetization curves at 2 K under 
cr, as shown in Fig. Efb). Although the changes in the 
M-H curve were rather small, we succeeded in observing 
the ferromagnetic hysteresis loops under cr by subtracting 
the data measured at cr = 1 MPa from the other data, 
as shown in Fig. He). The remanent magnetization and 
coercive field increases and decreases with increasing cr, 
respectively. 

Figure [31(a) shows the results of pyroelectric measure¬ 
ments. Before these measurements, the sample was 
cooled under applied cr, H = 3 T and a poling electric 
field ( E p ) of 222 kV/m, and then H and E p were removed 
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FIG. 2: (Color online) (a) M-H curves observed at 6 K under 
applied a of up to 151 MPa. The insets show (left-top) a 
schematic showing the directions of the Fe spins in the AFM 
phase, a and H , and (right-bottom) a magnification of the 
hysteresis loops near H — 0. (b) M-H curves and (c) their 
differences from the data at a — 1 MPa measured at 2 K 
under applied cr of up to 151 MPa. 


at 2 K. Subsequently, we measured pyroelectric current 
on heating under a. We found that the application of a 
induces spontaneous electric polarization along the c axis 
below T^ y , and that the magnitude of P c increases with 
increasing cr. We also observed that the sign of P c was 
reversed when the sample was cooled in a negative E p , 
as shown in the inset of Fig. [3fa). It is worth mention¬ 
ing here that in the previous study on Ba 2 CoGe 207 HU, 
polarity of the piezoelectric polarization was determined 
only by direction of the uniaxial stress. By contrast, in 
DyFeC> 3 , the application of cr induces a genuine ‘ferroelec¬ 
tric’ state, in which the spontaneous electric polarization 
can be reversed by an application of electric field. 

To confirm the that cr-induced ferroelectricity is of spin 
origin, we investigated H dependence of P c as follows. 
We cooled the sample under applied cr = 101 MPa, H = 3 
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FIG. 3: (Color online) (a) Temperature dependence of P c 
measured on heating in zero magnetic field after cooling under 
applied cr, H = 3 T and E v — 222 kV/m. The inset shows 
the data measured at a — 101 MPa after cooling with positive 
and negative poling electric fields, (c) and (e) show the values 
of P c at 2 K measured by pyroelectric measurements in zero 
magnetic field after temperature and magnetic field sequences 
shown in (b) and (d), respectively. 


T and E p = 222 kV/m. After removing E p at 2 K, we 
changed H from 3 T to a negative magnetic field of # 2 , 
and then removed it, as shown in Fig. [3fb). Figure [3fc) 
shows H 2 dependence of the values of P c at 2 K mea¬ 
sured by the pyroelectric measurements after the above 
mentioned sequences. We found that the cr-induced P c 
was flipped by the application of H 2 below ^ -2 T. We 
also performed similar magnetic field sweeping sequences 
shown in Fig. [3](d) , revealing that the flipped electric 
polarization was reversed again by applying a positive 
magnetic field of H% above ^ 2 T, as shown in Fig. [3](e) . 
These results show that P c shows hysteresis loop with re¬ 
spect to the magnetic fields applied at 2 K. Moreover, the 
width of the hysteresis loop roughly agrees with that of 
the M-H curves shown in Fig. EJc). Although the agree¬ 
ment between these hysteresis loops is not perfect, this 
is probably because we measured the remanent electric 
polarization at H = 0 after sweeping magnetic field up 
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FIG. 4: (Color online) (a) Piezoelectric responses at T = 2 K 
in zero magnetic field measured after cooling under a — 101 
MPa, H — 3 T and E p d= 222 kV/m. (b) Piezomagnetic 
responses at T = 2 K in zero magnetic field measured after 
H -sweeping from 0 to ±5 T and then to 0 T under a = 101 
MPa. 




to a certain value, while we measured the magnetization 
with varying H. From these results, we concluded that 
P c is coupled with the weak ferromagnetic moment, re¬ 
vealing that the cr-induced ferroelectricity originates from 
the spin order. 

Figure SJa) shows piezoelectric responses at 2 K in zero 
magnetic field. Before the measurements, the sample was 
cooled under a = 101 MPa, H = 3 T and E p ±222 kV/m. 
We removed H and at 2 K, and then measured the 
displacement current with varying a from 101 to 1 and 
then back to 101 MPa. The absolute values of P c were 
determined by the subsequent pyroelectric measurements 
on heating. We observed that the cr-induced P c changes 
in accordance with the varying cr, and that the polarity 
of the P c changes depending on E p applied on cooling. 

Similarly, we also confirmed the piezomagnetic nature 
as shown in Fig. |4jb). We swept H from 0 to 5 (or —5), 
and again to 0 T under cr, and subsequently measured 
the spontaneous magnetization at E[ = 0 with varying 
cr from 101 to 1, and again to 101 MPa. We observed 
that the cr-induced ferromagnetic moment also shows a 
similar behavior to that of the cr-induced P c . Note that 
the hystereses in these piezoelectric and piezomagnetic 
responses are probably due to the imperfect reduction of 
the stress; because of the friction between the zirconia 
pistons and anvils used in the pressure cell, the actual 
value of cr could not follow the nominal load value in 
particular in the cr-decreasing processes. 

Finally, we discuss the microscopic picture of the uni¬ 
axial stress effects on the magnetic orders. As mentioned 
above, the disappearance of the weak ferromagnetism at 
Xsr is related to the rotation of the G component of Fe 
spins The present experiments demonstrate that 


an application of cr induces weak ferromagnetic moments 
in the AFM and Dy-ordered AFM phases. From these 
results, we suggest that the application of cr affects in¬ 
plane magnetic anisotropy of the Fe spins, so that they 
are slightly rotated from their original position. The ro¬ 
tation of the Fe spins can yield a small G x component, 
which induces the spontaneous electric polarization in the 
Dy-ordered AFM phase via the magnetostriction mech¬ 
anism, in the same manner as in the field-induced mul- 
tiferroic phase 18, 0. By comparing the values of M 
and P c in the Dy-ordered AFM phase under a with those 
in the magnetic-field-induced multiferroic phase, the ro¬ 
tation angle is roughly estimated to be 1 ^ 2 degrees at 
cr = 151 MPa. 

In summary, we have investigated uniaxial-stress ef¬ 
fects on magnetic and dielectric properties of DyFeCE, 
which has a centrosymmetric crystal structure. By means 
of magnetization and pyroelectric measurements, we re¬ 
vealed that an application of cr induces both the weak fer¬ 
romagnetism and ferroelectricity in the Dy-ordered AFM 
phase. We also demonstrated that the magnetic-field- 
induced reversal of the weak ferromagnetic moment is 
accompanied by the flip of the electric polarization, re¬ 
vealing that the two cr-induced ferroic orders are coupled 
with each other. As for the microscopic mechanism of the 
uniaxial-stress effects, we suggest that the application of 
cr results in the change in magnetic anisotropy, which 
leads to a rotation of the Fe spins. This rotation yields 
the small G x component of the Fe spins, which is rele¬ 
vant to the emergence of the weak ferromagnetism and 
also induces ferroelectricity through the magnetostriction 
mechanism. This work provides a new guiding princi¬ 
ple for designing spin-driven multiferroicity. By apply¬ 
ing anisotropic stress on an antiferromagnet, one can re¬ 
move a number of symmetry operations in the system. 
This may lead to magnetically and electrically polarized 
states, in which the degree of the polarization can be 
tuned by the applied uniaxial stress. This guideline is 


also applicable for thin film samples [26|-[28j, in which 
structural mismatch between the film and substrate of¬ 
ten results in significant lattice deformation. By utilizing 
the epitaxial strain from the substrate, one may possibly 
impart additional ferroic orders to the film. 
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